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Preface  (to  Part  I) 


This  report,  composed  of  four  sections,  provides  a  summary  of 
work  accomplished  under  this  contract  in  four  diverse  areas  of 
optical  transmission  and  processing.  Briefly,  these  sections 
consist  of: 

I.  A  description  of  optic  of  wavelength  multiplexing.  This 
section  includes  as  an  attachment  a  report  on  the  subject 
presented  at  the  Conference  on  the  Physics  of  Fiber  Optics 
(Chicago,  April  1980). 

II.  A  summary  of  techniques  for  the  fabrication  of  charge  coupled 
devices  (CCD);  with  emphasis  on  the  selection  and  evaluation 
of  wafers,  oxidation  of  wafers  and  the  preparation  of  fibers 
on  the  oxidized  wafers. 

III.  A  short  summary  of  work  which  was  done  on  the  development  of  a 
device  to  produce  short  duration,  high  energy  optical  pulses. 

IV.  The  theoretical  development  of  optical  time  domain  reflecto- 
metry  applied  to  optical  waveguides  with  weak  static  perturba¬ 
tions.  This  section  also  includes  as  an  attachment  a  report 
on  the  subject  presented  at  the  Chicago  Conference. 

Other  work  for  this  contract  was  done  to  enhance  the  collection 
and  evaluation  of  data  on  optical  properties  of  various  solid 
state  materials.  This  work  included  both  the  development  of 
hardware  and  software  devices  for  this  purpose.  Part  II  of  this 
report  describes  the  hardware  and  software  which  was  developed  to 
interface  a  Nicolet  1080  computer  to  the  GPIB  (IEEE-488) 
specifically  for  use  with  a  video  camera  with  digital  output  (via 
thE  GPIB).  Other  work  in  this  area  has  been  previously  reported 
on  in  various  quarterly  status  reports  and  Technical  Memoranda. 
Briefly,  this  work  consisted  of: 

(1)  The  construction  of  a  data  encoder  for  use  with  various 
spectrometers  with  analog  output.  This  encoder  is  used  with 
a  Sykes  Comm-Stor  (Communications/Storage)  unit  to  create 
files  of  data  which  are  later  transferred  to  the  CDC6600  and 
analyzed  using  softare  created  for  a  Digilab  FTS-14  spectro¬ 
meter  (see  Final  Report  RADC-TR-78-61  for  details  of  the 
FTS-14  software). 

(2)  Revision  of  the  FTS-14  data  reduction  and  analysis  software 
for  use  with  various  signal  sources  including  the  above 
mentioned  encoder.  This  revised  software  is  used  with  output 
from  the  data  encoder,  a  punched  paper  tape  recorder  as  well 
as  the  FTS-14. 


(3)  Investigation  into  the  generation  of  correlation  masks  using 
the  AFGL  CRT  plotter.  This  included  the  generation  of 
resolution  test  strips  and  corrected  resolution  test  strips 
based  on  correction  coefficients  found  from  the  uncorrected 
strips.  In  addition,  assembly  language  programs  were 
developed  for  the  Nicolet  computer  to  test  the  feasibility  of 
using  this  device  to  generate  correlation  masks  in  real  time. 

(4)  A  set  of  software  tools  for  signal  processing  and  display, 
labeled  collectively  as  SIGPRODS,  was  modified  for  use  with 
parametric  plots.  This  is  specifically  used  with  software 
which  is  used  to  model  the  operation  of  a  class  of  tunable 
lasers. 
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SECTION  I 


Fiber  Optic  Wavelength  Multiplexers 


During  the  duration  of  this  contract  we  concentrated  our  efforts 
on  developing  new  concepts  for  fiber  optic  wavelength  multiplexers. 

In  particular  we  were  interested  in  devices  which  are  easy  to 
fabricate  and  thus  economical  and  which  display  low  insertion 
loss.  This  led  us  to  the  idea  that  one  could  devise  a  single 
element  multiplexer.  The  literature  on  multiplexers,  not  a  very 
extensive  one,  was  searched.  Basically,  multiplexers  utilize  two 
elements:  a  focussing-collimating  lens  and  a  dispersive  element. 

The  dispersive  element  can  take  the  form  of  a  prism  but  prefer¬ 
ably  a  diffraction  grating  to  keep  the  size  of  the  device  small. 

One  can  also  utilize  for  channel  separation  a  number  of  narrow 
band  interference  filters  arranged  in  a  structure  which  allows 
separating  the  various  channels  (optical  carriers)  at  different 
output  ports. 

The  concept  of  a  single-element  multiplexer  triggered  our 
interest  in  Fresnel  Zone  Plates  (FZP)  because  FZP's  act  as 
focussing  diffraction  gratings  and  have  high  dispersive  power. 

The  focal  length  of  a  FZP  is  inverse  proportional  to  wavelength. 

Let  the  focal  length  be  f0  at  a  wavelength  of  A0  at  an  arbitrary 
wavelength  x  •  we  have 

f(A)  =  fg-jf  (1) 


Taking  the  logarithmic  derivative  of  (1),  we  obtain 

!£  =  -  (2) 
f  * 

And  thus  the  axial  separation  between  two  adjacent  channels  is 

*f  -  £  J£1  (3, 


where  AA  is  the  interchannel  separation,  typically  20-30nm. 


On  the  guestion  of  low  insertion  loss  of  the  demultiplexer,  one 
must  be  careful  to  select  the  proper  type  of  FZP. 


A  "black-and-white"  FZP  consisting  of  alternative  opaque  and 
transparent  zones  has  a  maximum  efficiency  of  10%,  that  is,  it 
introduces  an  insertion  loss  of  10dB. 
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A  phase  reversal  FZP,  in  which  all  zones  are  transparent  but  they 
introduce  a  differential  phase  retardation  of  *  between  adjacent 
zones,  has  a  maximum  efficiency  of  40%,  i.e.,  an  insertion  loss 
of  6dB. 

One  can  do  better  by  utilizing  a  phase  Fresnel  lens  or  as  it  has 
been  recently  designated,  a  blazed  Fresnel  Zone  Plate  (BFZP). 

Such  a  device  is  capable  of  the  dispersion  given  by  equations 
( 1 ) — ( 3 )  but  introduces  negligible  insertion  loss  over  a  20% 
bandwidth. 

The  real  problem  that  remains  to  be  solved  is  to  find  a  good 
method  for  fabricating  a  BFZP.  Interest  in  BFZP  exists  also  in 
the  X-ray  field  where  it  can  be  used  as  an  efficient  lens.  By 
holographic  methods  one  can  create  FZP's  which  are  highly  effi¬ 
cient  in  materials  such  as  dichromated  gelatin.  The  problem  is 
that  most  of  the  interest  in  fiber  optics  is  in  the  near  infrared 
portion  of  the  spectrum  where  no  known  materials  exist  for  holo¬ 
graphy.  To  circumvent  this  problem,  one  can  create  the  FZP's 
with  an  argon  laser  at  a  visible  wavelength,  say  A  =  0.515um  and 
playback  in  the  IR.  Work  remains  to  be  done  to  reduce  the  aber¬ 
rations  which  occur  when  a  hologram  is  made  at  one  wavelength  and 
played  back  at  a  different  wavelength. 

A  paper  on  the  Fresnel  Zone  Plate  demultiplexers  by  Haim  Haskal 
has  been  presented  at  the  Conference  on  the  Physics  of  Fiber 
Optics  in  Chicago  on  April  28,  1980.  The  full  paper  is  scheduled 
to  appear  in  the  proceedings  of  the  conference. 
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SECTION  II 


Charge  Coupled  Device  Fabrication  Steps 


I.  INTRODUCTION 

During  the  contract  period  my  effort  has  related  larqely  to  the 
preparation  and  processing  of  silicon  wafers  to  he  used  in  the 
fabrication  of  experimental  charge  couDled  devices  (CCD)(U  For 
the  most  part  this  work  has  involved  day  to  day  contact  with 
members  of  ESE,  particularly  Mr.  Irvin  Berman.  Periodic  written 
reoorts  (quarterly)  and  special  oral  and  written  reports  have 
been  submitted  to  PML  and  the  Contract  Monitor. 

Major  emphasis  has  been  on  the  selection  and  evaluation  of  wa¬ 
fers,  the  oxidation  of  wafers,  and  the  preoaration  of  films  of 
polycrystalline  silicon  and  Silicon  Nitride  on  the  oxidized 
wafers. 

Along  with  the  day  to  day  laboratory  tVDe  activities,  literature 
surveys  and  evaluations,  somehwat  in  the  form  of  annotated  biblio¬ 
graphies,  have  been  made.  These  studies  include: 

a)  Limitation  imposed  on  CCD's  due  to  the  ouality  of 
available  silicon  wafers. 

b)  Problems  related  to  the  SiC>2/Si  interfaces. 

c)  Front  side  and  back  side  lettering  techniaues. 

d)  Expetaxial  growth  mechanisms  for  polycrystalline  silicons, 
silicon  nitride  (Si3N4)  and  low  temperature  silicon 
dioxide  (SiC>2)  layers. 

e)  New  technology  that  mav  impact  CCD. 

These  literature  surveys  (together  with  copies  of  almost  all  the 

articles  referred  to  therein)  are  available  in  ESE.  In  this 

report  we  are  merely  summarizing  some  of  the  conclusions  and/or 

recommendations  that  have  evolved  from  the  combination  of  the 

laboratory  activities  and  the  current  literature  studies.  These 

recommendations  and  conclusions  are  soecific  to  the  specialized 

ESE  environment  where  it  is  desired  to  make  sample  quantities  of  0 

research  devices  at  or  beyond  the  present  state  of  the  art.  The 

referneces  given  in  this  report  are  for  the  most  part  quite 

general  and  are  for  background  to  the  specific  area  rather  than 

justification  for  the  conclusions. 

* 
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II.  LIMITATIONS  IMPOSED  ON  CCD  DUE  TO  THE  QUALITY  OF  THE  WAFEP  USED 


1.  There  is  little  liklihood  that  any  material  will  replace 
silicon  as  the  semiconductor  from  which  CCD's  are  made.  Hess 
and  Sah (2)  show  that  the  hiqher  mobilities  of  less  tightly 
bonded  semiconductors  (e.q.,  III-V  Compounds)  cannot  be  effec¬ 
tively  used  because  of  the  saturation  of  the  drift  velocitv  of 
the  electrons  and  the  hot  electron  effects  in  the  fringing 
fields  inherent  in  the  CCD  geometry.  Jastrzebski  (3)  shows 
that  deep  recombination  centers  are  the  limiting  crystal  de¬ 
fect  in  silicon  CCD  devices.  These  types  of  defects  are  far 
more  numberous  in  compound  semiconductors.  Hence  this  reDort 
will  discuss  silicon  only. 

2.  Charge  Coupled  Devices  are:  extremely  sensitive  to  defect  intro¬ 
duced  both  in  the  silicon  growth  procedures  and  in  the  process¬ 
ing  steos  applied  to  the  silicon  wafers  prepared  from  the 
grown  crystals.  Hence  the  problems  for  achieving  the  objec¬ 
tives  of  Very  Large  Scale  Integration  (VLSI)  and  Very  High 
Speed  Integration  (VHSI)  are  closelv  related  to  the  corres¬ 
ponding  problems  in  CCDS.  Specifically,  these  objectives  are 
higher  density  of  smaller  elements  on  larger  chips,  higher 
speed,  greater  reliability,  hiqher  vields  and  lower  costs. 

3.  The  rapid  progress  made  in  the  last  few  years  on  the  objec¬ 
tives  cited  above  have  been  largely  due  to  improvements  in  the 
processing  steps  on  the  wafers.  There  has  been  considerable 
discussion  and  difference  of  opinion  on  whether  or  not  the 
present  quality  of  crystals  grown  by  the  Czochralski  or  Float 
Zone  technique  are  adequate  to  meet  the  needs  of  submicrometer 
size  devices.  Stated  another  way,  the  question  is  as  follows: 
can  improvements  (lower  temperatures)  in  the  processing  steos 
and  qettering  and  annealing  techniques  compensate  adequately 
for  the  defects  in  the  "as  grown"  crystals? 

First  we  will  briefly  summarize  tne  crystal  growth  situations. 
Practically  all  of  the  research  and  development  on  production 
and  experimental  devices  utilize  dislocation  free  Czochralski 
(CZ)  or  Float  Zone  (FZ)  crystals.  Consequently,  we  will  con¬ 
fine  our  attention  to  these  forms  of  silicon  crystal  growth. 

The  CZ  crystals  have  the  advantages  of  hiqher  physical  per¬ 
fection,  readv  availability  and  better  cost  and  yield  factors. 
The  FZ  crystals  have  greater  chemical  ourity  but  suffer  higher 
levels  of  physical  imperfections.  The  greater  chemical  purity 
of  the  FZ  crystals  is  due  to  the  elimination  of  the  quartz 
crucible  and  carbon  susceptors  and  resistance  heaters  which 
are  the  principal  source  of  oxvqen  and  carbon  impurities  in 
the  CZ  process.  The  greater  physical  perfection  of  the  CZ 
crystals  is  due  to  a  larger  melt  in  which  the  thermal  con¬ 
vection  currents  are  better  controlled.  Extensive  discussions 
of  the  status  of  CZ  and  FZ  growth,  the  effects  of  oxygen  and 
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carbon  and  the  types  and  nature  of  crystal  defects  are  given 
in  oublications  of  the  Electrochemical  Society and  the 
International  Conference  on  Crvstal  Growth  (ICCG)(*). 
excellent  review  of  the  limitations  imposed  on  the  CZ  growth 
by  thermal  convection  effects  has  been  published  by  J.R. 

Car r u ther s  )  and  by  Carruthers  and  Witt^).  A  similar  review 
emphasizing  FZ  crystals  has  been  published  bv  deKock^®).  These 
publications  and  a  large  number  of  current  paners  adeauatlev 
attest  to  the  difficulty  of  improving  either  of  these  types  of 
crystal  growths  (CZ  or  FZ).  A  summary  of  selected  DroDerties 
of  single  crystal  silicon  together  with  the  present  caoabili- 
ties  and  requirements  for  VLSI  is  given  in  Table  I  which  is 
taken  from  Bullish).  This  table  is  originally  from  a  report, 
as  yet  unpublished,  bv  a  select  committee  of  The  National 
Academy  of  Science (10).  The  table  clearly  shows  the  need  for 
improvements  of  a  factor  of  ten  or  better  in  many  cases.  Such 
improvements  are  unlikely  in  the  near  future  in  large  diameter 
crystals . 


TABLE  I 

Selected  Properties  of  Single-Crystal  Silicon 


Property 

Current  Capabilities 
Czochralski  Float  Zone 

Reauirements 
for  VLSI 

Resistivity  n-type  (P),.fl.*cm 

1-40 

1-300 

100-400 

Resistivity  n-tyoe  (Sb),/v*cm 
Resistivity  p-type  (B),-fl-*cm 

0.005-10 

0.005-50 

1-300 

0.0002-0.02 

100-400 

Resistivity  gradient  (four -probe)  % 

5-10 

20 

<1 

Resistivity  microsegregation,  % 

10-15 

20-50 

<1 

Minority  carrier  lifetime,  s 

30-300 

50-500 

300-1000 

Oxygen,  ppma 

5-15 

Not  detected 

Uniform  and 

Carbon,  ppma 

1-5 

0.1-1 

controlled 

0.1 

Heavy-metal  impurities,  ppba 

1 

0.01 

0.001 

The  conclusion  from  a  fairly  extensive  study  of  many  of  these 
factors  is  simply  that  ESE  should  continue  to  use  wafers  of  £  2 
inches  in  diameter  with  a  strong  preference  for  FZ  crystals.  The 
smaller  diameters  allow  much  greater  crystal  perfection  due  to 
much  better  control  of  convection  factors.  The  FZ  crystals, 
despite  larger  resistivity  striations  have  much  smaller  varia¬ 
tions  in  generation  lifetime.  This  is  in  turn  due  to  the  lower 
level  of  deep  recombination  centers  made  possible  by  reduction  of 
oxygen  and  carbon  in  FZ  crystals. 
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4.  Since  the  number  of  wafers  used  in  an  experimental  program  is 
quite  small,  and  the  research  time  spent  per  wafer  is  high, 
the  cost  per  wafer  is  relatively  insignificant.  Hence  one 
should  acquire  the  best  possible  pedigreed  wafers  to  the  most 
exacting  specifications. 

The  incoming  wafers  should  be  extensively  tested  on  a  sample 
basis  and  processed  ins  small  batches  of  the  order  of  ten  at  a 
time.  The  desirable  incoming  tests  are  a)  X-ray  topography, 
b)  Phase  Contrast  microscopy  on  "Secco"  etched  wafers  and  c) 
Capaci tance-voltage  measurements  (C-V).  All  three  of  these 
tests  are  time  consuming  and  the  latter  two  are  destructive  as 
well  (i.e.,  further  processing  steps  would  be  required  before 
reusing  the  wafer).  Nevertheless,  physical  tests  like  (a)  and 
(b)  are  essential,  as  well  as  electrical  tests  like  (c)  in 
order  to  correlate  electrical  properties  of  the  experimental 
devices  with  physical  and  chemical  properties  of  the  wafer. 

ESE  has  excellent  capability  for  all  three  of  these  tests. 

5.  what  is  really  needed  is  a  rapid,  non-destructive  imaging 
technique  that  would  put  a  high  resolution  image  of  the  defect 
pattern  of  a  wafer  onto  a  cathode  ray  tube.  The  best 
approach  (H)  »  (12)  to  this  type  of  device  is  to  use  laser 
scanning  beams  to  measure  surface  and  interface  properties  of 
wafers  before  and  after  various  process  steps.  DiStefano's 
device  has  been  developed  with  partial  support  of  government 
funds.  Such  a  piece  of  measuring  equipment  would  be  a  great 
asset  to  ESE. 

Bullish  has  written  an  excellent  review  of  methods  of 
characterizing  semiconductor  material.  Another  excellent 
recent  review  of  the  present  status  of  silicon  for  electronic 
devices  is  given  by  Bradshaw  and  Goor issen . 

It  is  strongly  recommended  that  Dr.  W.  Murray  Bullis  (National 
Bureau  of  Standards)  be  consulted  both  on  the  appropriate 
wafer  tests  and  on  acceptable  sources  of  wafers. 


III.  PROBLEMS  RELATING  TO  SILICON  OXIDATION 

Having  selected  and  prepared  the  wafer,  the  next  step  is  to  put 
onto  it  a  layer  of  the  order  of  0.1  to  1.0  micron  (1000  to 
10,000  A°)  of  silicon  dioxide  (Sj^).  The  current  accepted 
practice  is  to  use  a  dry  thermal  oxidation  at  a  temperature  of 
the  order  of  1000C°.  The  time  and  temperature  cycle  for  these 
oxidations  is  such  that  many  defects  native  to  the  wafer  are 
developed  and  many  more  may  be  introduced  in  this  processing 
step.  Furthermore,  the  properties  of  the  interface  between  the 
silicon  and  the  silicon  dioxide  layer  (Si/Sj^)  often  varied 
widely.  These  properties  (fixed  Oxide  Charges  Of  and  Nf)  are 
very  important  in  many  CCD  devices. 
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Because  of  inadequate  knowledge  of  the  oxidation  mechanism,  the 
physics  and  chemistry  of  the  Si/Si02  interface,  the  effects  of 
thermal  annealing  on  the  defect  structure  and  inconsistency  in 
definitions  and  nomenclature,  it  was  difficult  to  follow  and 
evaluate  much  of  the  copious  literature  in  this  area.  The  situ¬ 
ation  is  now  becoming  much  clearer.  The  nomenclature  and  defini¬ 
tion  problem  has  been  clarified  by  Deal  and  his  Commi ttee ( 14 ) . 

The  Si/Si02  interface  problems  are  reviewed  in  a  comprehensive 
way  in  the  Proceedings  of  the  International  Conference  on  the 
Physics  of  SiO^  and  Its  I nter f aces ( 15 ) .  The  understanding  of  the 
roles  of  chlorine(16)  and  hydrogen(l?)  in  frontside  gettering  has 
permitted  much  more  consistent  data.  Backside  getter ing ( 18)  is  a 
technique  wherein  purposely  introduced  imperfections  (disloca¬ 
tions)  in  the  inactive  portion  of  the  silicon  wafer  act  as  sinks 
to  remove  imperfections  and  imourities  from  the  critical  active 
device  regions  of  the  wafer.  This  techniaue  (backside  gettering) 
is  probably  the  largest  single  advance  in  device  processing 
introduced  in  the  last  several  years  in  terms  of  increased  yield 
and  quality  of  devices. 

The  combination  of  frontside  gettering  (using  chlorine),  backside 
gettering  (using  phosphorous  diffusion)  and  controlled  annealing 
should  allow  the  processing  of  excellent  Si02  layers  with 
superior  Si/Si02  interfaces  onto  excellent  quality  substrates. 

The  particular  implementation  of  these  techniques  depends  on  the 
circumstances  of  the  particular  group  involved.  The  particular 
recommendations  for  ESE  are  as  follows: 

a)  The  backside  gettering,  using  the  phosphorous  diffusion 
technique  of  Rozgonyi'^®)  is  preferred.  ESE  has  also  the 
capability  for  backside  gettering  by  ion  implantation  and/ 
or  mechanical  abrading  (sand  blasting). 

b)  The  frontside  gettering  technique  of  Hattori ( r ( 20)  which 
uses  a  N2/hC1/C>2  preoxidat  ion  ^anneal  is  preferred.  How¬ 
ever,  Tr ichlorethylene  (TCE)(21)  should  be  used  in  place  of 
HC1  as  it  is  less  corrosive  and  easier  to  handle.  One  sig¬ 
nificant  advantage  of  this  preoxidation  anneal  is  that  it 
is  done  "in  situ"  immediately  prior  to  the  oxidation  and 
acts  as  a  final  cleaning  step  for  the  wafer.  The  impor¬ 
tance  of  the  proper  combination  of  frontside  gettering, 
backside  gettering  and  slow  annealing  on  the  defect  struc¬ 
ture  of  Float  Zone  wafers  is  clearly  shown  by  Rohatgi  and 
Rai-Choudhury (22 ) .  Their  technique,  which  can  be  readily 
combined  with  that  of  Hattori  yields  millisecond  carrier 
lifetime  measurements  which  are  indicative  of  very  high 
quality  defect  free  wafers.  These  combined  technioues 
applied  to  two  inch  diameter  Float  Zoned  wafers  should 
provide  high  quality  oxide  layers  with  good  interface 
properties  and  potentially  excellent  CCD  device  properties 
(low  leakage,  etc.). 
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IV.  DEPOSITION  OF  POLYSILICON  LAYERS 

The  preparation  of  thin  films  of  semiconductors  onto  substrates 
of  identical  composition  (epitaxy)  and  onto  substrates  of  differ¬ 
ent  composition  (hetero-epitaxy)  by  Chemical  Vapor  Deposition 
(CVD )  techniques  has  increased  very  rapidly  over  the  past  decade 
and  it  can  easily  be  predicted  that  this  trend  will  continue  for 
many  years  to  come.  There  are  several  reasons  for  this  situa- 
t  ion: 

a)  The  planar  technology  has  required  relatively  thin  active 
layers  with  precisely  controlled  properties.  The  drives 
toward  VLSI  and  VHSI  sharply  accentuate  the  importance  of 
thin  high  quality  layers  (usually  in  the  region  of  100  to 
1000  A°) . 

b)  CVD  is  an  extremely  versatile  technique  with  a  very  wide 
range  of  applications  for  elemental  and  compound  semi¬ 
conductors,  conductors,  resistors,  and  capacitative 
components. 

c)  The  cost  and  complexity  of  setting  up  an  epitaxial  system 
is  far  less  than  to  set  up  a  Czochralski  or  Float  zone 
facility. 

d)  There  are  very  many  opportunities  for  practical  applica¬ 
tion  investigations  as  well  as  fundamental  investigations 
of  CVD  processes.  In  other  words,  the  potential  field  for 
new  research  opportunities  is  much  larger  and  the  state  of 
the  art  is  far  less  developed  than  that  of  growth  from  the 
melt. 

e)  There  are  very  many  investigations  that  are  suitable  for 
doctoral  theses  in  universities  without  excessive  invest¬ 
ment  in  capital  equipment.  The  individuals  thus  trained 
are  highly  desirable  in  small  and  specialized,  as  well  as 
large  and  generalized  industrial  organizations. 

f)  New  epitaxial  techniques  will  probably  replace  many 
presently  used  processing  steps  (see  Section  V  of  this 
report  -  New  Technologies). 

There  are  a  number  of  excellent  reviews  on  the  subject  of 
Chemical  Vapor  Deposition  of  semiconductors  (23)  and  particularly 
of  silicon(24)<  Qur  attention  has  been  confined  to  the  depo¬ 
sition  of  polysilicon  by  CVD  with  emphasis  on  low  pressure  CVD 
(LPCVD)  and  plasma  enhanced  CVD  (PECVD). 

The  advantages  of  LPCVD  are  clearly  spelled  out  by  Hammond(23)# 

He  outlines  a  simplified  barrier  layer  theory  in  which  he  shows 
that  the  diffusion  rates  for  the  gaseous  species  are  inversely 
proportional  to  the  total  system  pressure  and  directly  propor¬ 
tional  to  the  partial  pressure  of  the  gaseous  species  in  the 
mainstream.  Also  the  diffusion  rate  is  directly  proportional  to 
temperature  to  the  three-halves  power  (tV2).  Thus,  by  reducing 
pressure  from  one  atmosphere  to  a  fraction  of  a  TORR,  the 
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diffusion  rate  is  increased  by  a  factor  of  about  a  thousand  while 
the  deposition  rate  is  reduced  very  little  since  no  dilutent 
gases  are  used.  Thus,  the  uniformity  of  the  deposited  layers  is 
greatly  enhanced,  within  the  individual  wafers,  from  wafer  to 
wafer,  and  from  batch  to  batch.  The  merits  and  details  of  LPCVD 
are  discussed  by  Rossler^4)f  by  ogirima  et  al  ( 25)  r  and  by  Spear 
and  Wang ( 26 ) # 

The  specific  applications  of  LPCVD  for  ESE  arc  discussed  in  an 
internal  report (27) #  This  report  discusses  two  LPCVD  systems. 

One  of  these  is  essentially  a  "turn  key"  retrofit  system  to  be 
installed  on  the  ESE  Lindberg  furnace.  The  other  is  an  "in 
house"  experimental  system  which  is  scaled  down  for  specific 
needs.  Operation  of  the  first  system  is  now  scheduled  for  April 
1981.  The  second  system  is  currently  being  "debugged". 

It  is  recognized  that  Plasma  Enhanced  CVD  systems  (PECVD)  may 
have  many  advantages  over  LPCVD  systems.  However,  they  also  have 
additional  compilations  so  that  there  will  be  little  emphasis  in 
ESE  on  PECVD  until  after  the  LPCVD  is  fully  implemented.  It  is 
the  present  intent  to  use  a  system  which  is  a  modification  of 
that  Rosier  and  Engle(28)4  The  current  thinking  concerning  the 
proposed  system  is  discussed  in  the  internal  reoort  cited (27) # 

It  represents  a  conversion  of  the  present  MiniPak  II  system  to  an 
R.F.  Heated  Hot  wall  system. 


V.  NEW  TECHNOLOGIES  THAT  MAY  INFLUENCE  NEW  CCD  DEVICES 

Some  effort  has  been  expended  in  exploring  new  technologies  that 
may  influence  future  directions  of  VLSE,  VHSI  and  CCD  research. 
They  will  merely  be  mentioned  here.  The  most  significant  trend 
in  my  opinion  is  the  increasing  significance  of  Chemical  Vapor 
Deposition.  A  good  part  of  the  recent  success  has  related  to 
better  control  of  the  variables  involved  and  to  higher  purity 
intermediate  starting  materials  such  as  SiCl4,  SiHCl3  and  SiH4« 

The  major  limitations  of  CVD  has  been  the  substrates.  The  possi¬ 
bilities  of  reusing  substrates  by  cleaving  the  epilavers,  the  use 
of  contoured  foreign  substrates  (g raphoepitary )  and  the  use  of 
laser,  electrons  and  ion  beam  annealing  of  epitaxial  layer  or 
foreign  substrates  may  radically  change  our  current  processing 
technologies  in  the  next  decade.  The  ability  to  deposit  Si02 
layers  by  CVD  techniques  whose  properties  are  comparable  or 
superior  to  thermal  Si02  is  apparently  possible < 2^) . 


Despite  theoretical  arguments  that  ideal  Si02-Si  interfaces  are 
possible ( 30) t  it  is  unlikely,  in  my  opinion,  that  the  fixed 
charge  density  at  the  Si02/Si  interface  will  be  reduced  much 
below  the  present  best  values  of  the  order  of  1010cm~2.  Despite 
the  central  importance  that  Si02  has  had  on  the  development  of 


our  current  electronic  devices  and  systems,  it  may  be  necessary 
to  develop  a  new  insulating  interface  with  appropriate  crystal¬ 
lographic,  as  well  as  thermal  and  mechanical  properties  if  we  are 
to  greatly  improve  on  present  interface  conditions  and  to  make 
truly  three-dimensional  (layered)  VLSI  and  VHSI  devices.  One 
possibility  would  be  magnesium  Spinel  MgAl204  if  it  could  be 
deposited  on  silicon  substrates  in  single  crystal  form  by 
metaloorganic  techniques.  Then  multiple  alternate  layers  of 
single  crystal  silicon  and  spinel  could  be  used  to  explore  the 
feasibility  and  desirability  of  truly  three-dimensional  devices. 

The  rapid  development  of  computer  aided  design  and  computer  aided 
manufacturing  techniques  (CAD/CAM)  will  greatly  accelerate  new 
VLSE,  VHSI  and  CCD  design  and  prototype  chips  at  greatly  reduced 
cost . 


VI .  FUTURE  DIRECTIONS 

It  is  anticipated  that  next  year's  efforts  will  be  much  more 
restricted  in  scope  and  relate  the  experimental  aspects  of  the 
LPCVD  and  PECVD  programs. 
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SECTION  III 


Generation  of  Short,  High-Energy  Pulses 

During  this  period,  a  major  effort  was  directed  toward  develop¬ 
ment  of  a  circuit  to  produce  repetitive,  short,  high-energy 
optical  pulses.  The  initial  design  goal  was  an  optical  pulse 
duration  of  about  10  ns.  Given  sufficient  energy  per  pulse,  an 
input  pulse  of  this  duration  will  significantly  enhance  the 
resolution  performance  of  the  Optical  Time  Domain  Pef lectometer 
presently  in  operation. 

The  circuit  used  in  this  developmental  optical  pulser  employs  a 
charged  delay  line  as  the  energy  storage  element,  with  an  n-p-n 
transistor  operating  in  the  avalanche  mode  as  a  switch  connecting 
the  delay  line  to  a  GaAs  laser  diode.  A  series  of  experiments 
were  performed,  using  both  PC*  &G2000  series  laser  diodes  and  the 
Laser  Diode  Laboratories  .eries  diodes.  The  choice  of  the 

switching  transistor  turfis  u not  to  be  critical,  so  long  as  the 
reverse  collector-to-ba*»  >.«.'< -  vaown  voltage  is  greater  than  the 
400  to  500  volts  needed  i-o  *  *  ge  the  delay  line.  With  this 

configuration,  using  a  M'  t'.o^a  2N5550  transistor  and  microstrip 
delay  lines  of  characve*  is’:  j  ;  impedances  between  2  and  10  ohms, 

15  ns  optical  pulses  *»£  up  lo  20  watts  Deak  power  at  900  nm  can 
be  reliably  obtained  at  rates  of  1000  pps. 

The  experiments  described  above  employed  discrete  circuit  fabri¬ 
cation  techniques.  Using  these  techniaues,  it  is  possible  that 
optical  pulses  as  short  as  5  ns  duration  could  be  generated  at 
the  desired  power  levels.  The  limitation  on  the  use  of  discrete 
circuit  components  results  from  the  effects  of  distributed  induc¬ 
tance  throughout  the  circuit.  To  obtain  a  feel  for  the  magnitude 
of  these  effects,  at  a  frequency  of  170  MHz,  a  1-centimeter 
length  of  conductor  has  a  reactance  of  about  10  ohms;  to  produce 
20  watt  optical  pulses  using  existing  GaAs  laser  diodes,  the 
driving  pulse  current  approaches  100  amperes.  At  these  current 
levels,  the  inductive  voltage  drop  across  even  a  few  millimeters 
of  conductor  becomes  substantial  as  the  rise  time  of  the  current 
driving  pulse  approaches  a  few  nanoseconds. 

The  most  promising  next  step  appears  to  be  implementation  of  this 
circuit  in  hybrid  form.  C.  Luddington,  of  ESED  in  the  Solid 
State  Sciences  Division,  fabricated  0.5  ohm  resistors  on  0.01 
ohm-centimeter,  single-crystal  silicon  wafer  material  and  mounted 
several  of  these  resistors  in  headers.  Measurements  of  these 
encourage  the  notion  that  a  hybrid  circuit,  with  the  switch 
transistor,  current-limiting  resistor  and  the  laser  diode  all 
mounted  in  the  same  header,  would  greatly  reduce  distributed 
inductance  in  the  circuit,  hence  making  possible  shorter  optical 
pulses. 
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SECTION  IV 


Theory  of  Optical  Time  Domain  Ref lectometry 
From  Perturbed  Fiber's 


We  have  developed  a  theory  to  describe  the  Optical  Time  Domain 
Ref lectometry  (OTDR)  return  signals  from  weak  static  perturba¬ 
tions  in  optical  waveguides  by  employing  a  time-dependent  Green's 
function  formulation.  The  version  of  this  theory  that  is  appli¬ 
cable  to  two-dimensional  Dlanar  waveguides  has  been  presented  at 
the  American  Ceramic  Society's  1980  Chicago  meeting  and  will 
appear  in  Advances  II:  Physics  of  Fiber  Optics.  It  is  included 
as  part  of  this  report^!. 


Further  calculations  for  the  planar  waveguide  case  have  been 
carried  out  in  order  to  establish  qualitative  comparison  with 
experiments  that  have  been  performed  for  perturbations  consisting 
of  loops  12)  of  various  diameters  ranging  from  5  to  12  cm.  To 
facilitate  the  numerical  evaluation  of  the  theory,  we  have  car¬ 
ried  out  the  integration  of  Equation  37  of  Ref.l  analytically. 
This  yields  an  expression  for  the  backscattered  Rayleigh  oower, 
that  can  be  written  as: 
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where : 


£  '  2  (  to  Tot  ~  *  Vj~'  /  ZX  ) 

]  '  2(«tox  -tvj'/z1*  oc.gAt^ 


In  Figure  1.  we  have  plotted  the  local  attenuation  in  dB/Km  vs. 
the  actual  position  along  the  waveguide  for  loops  of  various 
diameters  for  the  parameters  of  the  planar  model  of  Ref.l  using 
the  expression  for  the  Rayleigh  backscatter ing  power  given  above 
in  Equation  1. 


As  mentioned  on  page  7  of  Ref.l,  inhomogeneities  on  a  microscopic 
scale  are  essentially  the  mechanism  that  produces  Rayleigh  scat¬ 
tering  and,  for  the  case  of  glass  fibers,  these  inhomogeneities 
are  due  to  thermal  fluctuations  in  density  and  composition  that 
become  frozen  in  as  the  glass  melt  solidifies  at  temperature  Tc. 
These  local  fluctuations  of  density  and  composition  will  in  turn 
result  in  local  fluctuations  of  the  refractive  index.  A  calcu¬ 
lation  of  Rayleigh  scattering  can  thus  be  carried  out  using  the 
Green's  function  methods  discussed  in  Ref.l,  with  the  perturba¬ 
tion  in  refractive  index  described  by 


where  A ya(r)  is  the  fluctuation  in  density  at  point  r  . 
The  scattered  field  can  then  be  expressed  as 
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where  |w>  is  the  matrix  element  of  the  density  fluctua¬ 


tion 


oo 


<°  l  ^ o.  1  n  >  *  f  uo O )  zy)  (  *,  u  „  i x) 

’  —  oo 


y ) 


between  the  incoming  mode  U.  and  the  scattered  mode  U.„  ,  and  i 
labels  the  fluctuation  volume  which  is  centered  around  Zj  where 


zi+l  ~  zi  ~  ^ 


The  total  back-scattered  power  relative  to  the  input  oower 
will  be  given  by 
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where 
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«j  tot 

and  the  bar  indicates  a  thermal  average. 
Since 


and 


(7) 


the  backscattered  power  from  the  extrinsic  perturbation  and  the 
intrinsic  frozen  in  thermal  fluctuations  will  not  interfere  and 
the  total  backscattered  power  will  be  simply  the  sum  of  the  two 
component  powers.  The  relative  backscattered  Rayleigh  oower  can 
be  written  as 

(6) 
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Instead  of  pursuing  a  detailed  evaluation  of  this 
we  shall  obtain  an  approximate  value  for  it  in  a 
offer  a  comparison  with  the  formula  derived  using 
To  do  this  we  first  approximate  all  ^n(and  >7*)  as 
the  s_.ne  and  equal  to  /3C  (and  -u,  )  .  Having  made  th 
tion,  we  next  assume  that  enough  guided  modes  exi 
sum  over  intermediate  states  can  be  approximated 
operator.  The  quantity  a  can  then  be 

thermodynamics  1*1  to  be  /  /  ' 


quantity  here, 
rough  way  to 
ray  theory  ^  . 
being  roughly 
is  approxima- 
st  so  that  the 
by  the  identity 
calculated  by 


(9) 


^  =  kT/3r  yOVV 


which  we  further  assume  to  be  roughly  constant  over  the  guide 
volume.  For  purposes  of  comparison,  we  take  the  guide  to  be 
planar  so  that  the  fluctuation  volume  can  be  roughly  taken  as 
V/*  A#  ha*  with  a  the  guide  width.  For  a  multimode  step  index 
guide,  we  approximate  R.Vas 
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which  then  yields  a  rough  upper  bound  for  Rayleigh  backscatter ing 
in  the  guide 
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The  Green's  function  formulation  for  OTDR  in  planar  waveguides 
may  be  extended  to  the  case  of  circular  fiber  guides  in  a 
straightforward  manner.^  Replacing  the  eigenfunctions  1L;  of 
the  planar  guide  by  \Aj  yfr)  where  l,n  designate  the  angular  and 
radial  mode  numbers  with  r  the  radial  coordinate,  the  Green's 
function  for  the  circular  case  may  be  written  as 
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(13) 


where  the  summation  symbol  indicates  an  integral  continuum  state 
as  well  as  summation  over  the  bound  modes  which  we  abbreviate  as 
•  For  &n  initial  pulse  with  a  modal  content  described 
by  the  vector  aj,  i.e.. 
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scattering  in  the  backward  or  forward  direction  into  the  k ^  mode 
will  be  given  by 

"  Kw  e' '  fa  )  a  -  $*£  n]  rjz 


J5) 


where 
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=  2ir/  u*  [ShV)]  U-c**)  rdv 

and  where  we  have  assumed  an  abrupt  transition  for  the  oerturbed 
region  starting  at  z  =  a  and  extending  to  z  =  a+L  .  For  this 
case 
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The  expressions  for  transition  loss  and  reflected  power  are 
sentially  the  same  as  those  for  the  planar  case  with  the  set 
=  replacing  the  transverse  mode  number  n  in  the  planar 

case.  As  an  example,  the  expression  for  the  power  lost  due  to 
transition  radiation  will  be  given  by 
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The  transverse  eigenfunctions  for  the  steo  index  circular  guide 
are  more  complex  than  those  for  the  planar  guide,  with  Bessel 
functions  replacing  simple  exponentials^.  For  the  case  of  the 
bend  perturbation,  the  double  sum  in  the  expressions  for 
scattered  power  can  be  simplified  due  to  the  selection  rule 
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which  allows  transitions  only  to  final  states  A  that  differ  from 
the  initial  state  J  by  one  unit  of  the  angular  mode  number. 


The  radial  portion  of  the  transverse  matrix  element  may  be 
evaluated  by  reexpressing  the  integrals  involving  products  of 
Bessel  functions  as 
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with  similar  expressions  for  the  products  of  the  other  Bessel 
functions.  The  complete  list  of  integrals  needed  for  the 
calculation  of  for  the  case  of  bends  is  9iven  in  Appendix  A. 
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The  only  other  remaining  difference  between  the  circular  and 
planar  calculation  for  OTDR  from  bends  comes  from  the  expression 
for  radiation  loss  due  to  tunneling.  For  the  circular  case,  this 
is  given  by^6^ 
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which  has  the  same  exponential  dependence  upon  the  bend  radius  R 
as  does  the  two-dimensional  case  plus  an  additional  factor  r!/2 
in  the  denominator. 

A  computer  code  to  evaluate  OTDR  returns  from  bends  in  circular 
fibers  has  been  prepared  and  is  in  the  process  of  being  debugged. 
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APPENDIX  A 


The  integrals  needed  to  evaluate  the  transverse  matrix  element 
for  the  case  of  bends  in  a  step  index  fiber  are  tabulated  below 
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Using  these  results,  the  transverse  matrix  element  for  bound  to 
continuum  transitions  may  be  expressed  (writing  *  for  r**,  <■  and 
7^  for  y>  ;  )  as 
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ATTACHMENT  A 


FRESNEL  ZONE  PLATES  FOR  WAVELENGTH  DEMULTIPLEXING 


by 

Haim  M.  Haskal* 

Parke  Mathematical  Laboratories 
Carlisle,  MA.  01741 


An  analysis  is  presented  on  the  application  of  Fresnel  zone 
plates  to  wavelength  demultiplexing  in  fiber  optic  systems.  A  de¬ 
multiplexer  consisting  of  an  off-axis  Fresnel  zone  plate  acts  as  a 
focusing  spectrometer.  Blazed  Fresnel  zone  plates  have  high  efficiency 
over  a  20‘X,  bandwidth  and  can  demultiplex  a  large  number  of  channels. 


* Research  supported  by  Solid  State  Sciences  Division,  Rome  Air  Development 
Center  (AFSC)  under  contract  No.  F19628-78-C-0089. 
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INTRODUCTION 

Present  day  optical  fibers  have  capacities  which  exceed  by  far  the 
bandwidths  required  to  transmit  high  speed  modulated  light  signals.  To 
increase  the  bandwidth  utilization  of  the  fiber  several  light  carriers  can 
be  simultaneously  transmitted  -  multiplexed  -  over  the  same  fiber  and  separ¬ 
ated  by  wavelength  -  demultiplexed  -  at  the  receiving  end.  The  first  propo¬ 
nent  of  wavelength  multiplexing  was  Tomlinson  at  Bell  Telephone  Laboratories 
who  ab.lyzed  ^  in  detail  the  requirements  and  techniques  for  optical 
multiplexing.  Optical  demultiplexers  based  on  blazed  diffraction  gratings 
were  demonstrated  by  Tomlinson,^ and  Aoyama^.  Other  types  of  demulti¬ 
plexers  based  on  interference  filters^,  prisms^,  and  aberration  corrected 
concave  gratings^  were  also  demonstrated. 

In  this  paper  we  analyze  the  performance  of  a  wavelength  demultiplexer 
based  on  a  Fresnel  zone  plate.  The  zone  plate  performs  simultaneously  the 
functionsof  wavelength  dispersion  and  focusing. 

FRESNEL  ZONE  PLATE  THEORY 

Fresnel  zone  plates  (FZP)  have  been  used  as  lenses  and  as  spectro¬ 
meters  in  various  regions  of  the  electromagnetic:  spectrum.  An  excellent 
article  reviewing  the  properties  of  FZP  was  published  by  Young  The 
usual  type  FZP  consists  of  a  series  of  alternating  opaque  and  transparent 
annular  zones.  This  FZP  is  characterized  by  a  transmission  function  (for  an 
infinite  plate) 


t(r) 


1 
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sgn  (cos 


kr2 

FT5 


0) 


where  r  is  the  radial  coordinate  in  the  FZP  plane,  k  =  2ir/\,  \  1S  the 
wavelength,  f  the  focal  length. 

For  a  normally  incident  unit  plane  wave  the  transmitted  electric 
field  can  be  expanded  in  a  Fourier  series 

“  sin  not  .  2 

E(z,r)  =  exp(-jkz)  L  T~  exp(j  --) 

n=-«>  nit 
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In  Eq.  (2)  the  terms  with  n_ negative  represent  spherical  waves  con¬ 
verging  toward  real  foci  of  coordinates{0, ^ ) .  Similarly  the  terms  with  n 

positive  represents  spherical  waves  diverging  from  virtual  foci  located  at 


coordinates  (0,  "-j-jj-  )  as  shown  in  Fig.  1.  The  fraction  of  the  incident  light 

diverted  to  the  first  real  focus  is  called  the  efficiency  of  the  FZP  and  in 
this  case  is  given  by 

n  =  -5  =  10%. 

*]  IT  * 


A  more  efficient  zone  plate  can  be  constructed  by  replacing  the  opaque 
zones,  which  obstruct  50%  of  the  Incident  light,  by  transparent  zones 
which  introduce  a  half  wave  retardation  between  adjacent  zones  as  shown  In 
Fig.  2.  Such  a  zone  plate  is  called  a  phase  reversal  zone  plate  and  is  also 
characterized  by  multiple  foci  but  has  an  efficiency  n  given  by 
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2  40  X. 


The  focal  length  f_  for  both  types  of  Fresnel  zone  plates  mentioned 

t  h 

above  can  be  related  to  the  zone  number  m  and  the  radius  of  the  m  zone  by 


f  = 


mX 


(3) 


Let  the  focal  length  of  the  FZP  be  f  at  a  wavelength  XQ  then  at  an  arbitrary 
wavelength  _X_ we  have 

f(A)  *  f„  f  (4) 

Equation  (4)  puts  in  evidence  the  high  dispersion  of  FZP's  which  can  be  used 
advantageosuly  in  wavelength  demultiplexing. 

To  improve  the  efficiency  of  the  FZP  beyond  the  40%  that  is  possible 

with  the  phase  reversal  zone  plate  one  can  use  the  phase  Fresnel  lens  concept 

(91 

introduced  by  Miyamoto'  '  in  1961.  Instead  of  introducing  an  abrupt  phase 
retardation  of  a/2  at  the  edge  of  a  zone  as  in  the  phase  reversal  zone  plate, 
in  Miyamoto's  phase  Fresnel  lens  one  introduces  a  gradual  phase  retardation 
over  a  zone,  according  to  some  profile,  from  zero  a  maximum  of  X.  The 
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adjacent  zone  starts  again  from  a  zero  phase  retardation  which  Is  indistinguishable 
from  the  _A_  retardation  at  the  edge  of  the  previous  zone.  In  this  way  the  phase 
retardation  appears  continuous  for  the  incident  light  even  though  the  actual 
phase  plate  has  a  dentelated  profile.  The  procedure  of  constructing  a 
phase  Fresnel  lens  corresponding  to  a  plano-convex  lens  is  illustrated  in 
Fig.  3. 

The  phase  Fresnel  lens,  also  referred  to  as  the  blazed  Fresr.el  zone 
plate  (BFZP),as  designed  in  Fig.  3  will  focus  a  plane  wavefront  at  the  design 
wavelength  XQ  to  a  single  focus  just  like  the  original  lens.  However,  if 
the  wavelength  is  changed  from  \Q  to_x_some  of  the  light  is  focused  to  higher 
order  foci.  We  will  now  examine  the  effect  of  changing  the  wavelength  on  the 
efficiency  of  the  BFZP. 

Let  the  BFZP  introduce  a  quadratic  wavefront  deformation  given  by 
the  phase  function 
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Substituting  for  the  coefficients  in  eq.  (7)  we  obtain 


sin  u(x  -X)/X 

exp(- j2ti<j>/x)  ■  - - -  exp  {-j2no/X-jir( XQ-X) /X  } 

n(XQ-X)/X 


x  1  +  -2 —  exp(j2ir0/x)  ♦ 

L  x° 


exp(-j27re/x)  +  . 


J 


(8) 


At  the  design  wavelength  \Q  only  the  first  term  in  the  series  survives 
corresponding  to  a  singleTeal  focus.  As  one  deviates  from  XQ  higher  order 
foci  appear.  To  estimate  the  fraction  of  energy  diffracted  7t  the  first  real 
focus,  i.e.  the  efficiency  n,  as  a  function  of  wavelength,  we  consider  the  co 
efficient  of  the  first  term  in  Eq.(8).  The  efficiency  is  given  by 


'sin  2^(1  ♦  &)  1 
Ao  Ao 

'1 =  - 

(l  +M.) 

-  ^  o  *0  * 

where  AX  is  the  deviation  from  the  design  wavelength  XQ.  If  we  set  AX/ ^ 


(9) 


+_  10%  in  Eq.  (9)we  obtain  the  value  0.96;  thus  only  4%  of  the  energy  is 
diffracted  into  the  zero  and  higher  order  foci.  A  BFZP  can  therefore  be  an 
efficient  demultiplexer  with  20%  bandwidth.  Such  a  bandwidth  can  accommodate 
10  optical  channels  20  nm  apart  around  a  center  wavelength  ^  =  1.0  pm. 
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THE  USE  OF  THE  BFZP  FOR  WAVELENGTH  DEMULTIPLEXING 


The  BFZP  can  be  used  in  several  configurations  (reflection,  transmission) 
as  a  dispersive  and  focusing  device  for  wavelength  demultiplexing.  It  has, 
in  fact,  been  proposed  for  use  as  a  focusing  spectrometer  at  soft  x-ray 
wavelengths^).  One  cannot  use  a  full  BFZP  for  demultiplexing  because  the 
images  of  the  input  fiber  at  the  various  wavelengths  appear  along  the  axis 
one  in  the  shadow  of  the  other.  If  an  array  of  detectors  is  placed  at  the 
image  positions,  the  first  detector  will  obstruct  some  of  the  light  meant  to 
reach  the  others,  in  addition  crosstalk  between  channels  will  be 
very  high.  To  alleviate  this  problem  one  can  use  an  off-axis  BFZP;  light 

is  now  focused  towards  the  various  images  only  from  one  side  of  the  plate  and 
the  array  of  channel  detectors  can  be  placed  along  the  BFZP  axis  or  at  an 
angle  to  the  axis  without  interference. 

We  illustrate  the  use  of  a  BFZP  in  an  off-axis  configuration  in  Fig.  4. 
The  light  emanating  from  a  multimode  optical  fiber  is  collimated  prior  to 
entering  the  BFZP.  Two  adjacent  optical  channels  at  and  An  +  A A  are 
focused  at  F^  and  F?.  The  separation  between  the  two  channels  along  the 
axis  can  be  obtained  from  Eq.  4  and  is  given  by 


tf 


(10) 


where  Axis  the  wavelength  separation  between  the  channels.  From  tq.  (10)  one 
determines  the  size  of  photodetectors  as  a  function  of  the  BFZP  focal 
length  and  wavelength.  The  light  impinges  on  the  detectors  located  along 
the  axis  at  a  very  shallow  angle,  close  to  grazing  incidence  and  therefore  • 
the  reflection  loss  can  be  considerable.  It  is,  therefore,  necessary  to 
coat  the  photodetectors  with  antireflection  coating* suitable  for  that 
angle  of  incidence. 

We  have  shown  in  Fig.  5  a  two-channel  demultiplexer  with  the  detectors 
inclined  at  an  angle  6  to  the  BFZP  axis;  they  are  located  such  as  to  inter¬ 
cept  all  the  light  from  the  respective  image  of  the  fiber. 
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The  angle  6  is  determined  from  the  geometry 


tan  6 


r  -  D  +  ea 
m 

f 


(id 


where  rmis  the  radius  of  the  mth  zone  of  the  BFZP,  £  is  the  diameter  of  the 
collimaled  beam,  £  is  the  angle  of  divergence  of  the  collimated  beam,  £  is 
the  distance  between  the  collimating  lens  and  the  BFZP. 

We  also  have  the  following  relations 

0  =  2(NA)  fC£)1  (12) 

where  (NA)  is  the  numerical  aperture  of  the  fiber,  f^  the  focal  length 
of  the  collimating  lens. 


c 


core 


(13) 


where  dcQre  is  the  core  diameter  of  the  fiber. 

From  Eq.(3)we  obtain  the  relation 
1_  =  AH 

rm  A 


(14) 


where  r  is  the  width  of  the  mth  zone.  Substituting  into  Eq.  (11)  we  obtain* 


tan  <S 


A 

W 


-2(NA) 


(15) 
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where  we  have  assumed  that  the  term  ea  In  Eq.  (11  )is  negligible. 

i  col 

If  one  uses  a  unity  magnification  imaging  system  =  1;  NA  =  0.2,  X  -  )  um 

and  Ar  =  1  pm  we  obtain 

tan  <5  =  0.6  6  =  30.96° 


Thus  in  this  case  the  anti  reflection  coating  for  the  photodetector  must  be 
designed  for  light  incident  at  59°  with  respect  to  the  normal  if  the  detector  is 
placed  along  the  axis  of  the  BFZP.  If  the  detectors  are  inclined  at  angle 
6  with  respect  to  the  axis  the  antireflection  coating  should  be  designed 
for  an  incidence  angle  of  28°  with  respect  to  the  normal  to  the  photodetector. 

The  dimensions  of  the  BfZP  need  not  be  very  large  and  are  compatible  with 
fiber  optic  systems.  If  we  choose  f  =  fcoj  =  6  mm  ,A  =  1  pm  ,  Ar  =  1  pm,  the 

diameter  of  the  collimated  beam  equal  to  the  off-axis  BFZP  diameter  is 

0=2  (NA)fCo]  =  2.4  mm 

The  number  of  zones  which  must  be  ruled  is  readily  computed  to  be  1920. 

FABRICATION  CONSIDERATIONS 

To  fabricate  a  blazed  Fresnel  zone  plate  several  avenues  are  open: 

a)  Diamond  point  turning.  Considerable  progress  has  been  reported  in  the 
last  few  years  in  the  design  and  implementation  of  computer  controlled  diamond 
turning  machines.  Such  machines  are  capable  of  generating  the  required  phase 
profile  in  soft  plastic  materials  to  produce  a  BFZP. 

m 

b)  Computer  generated  phase  profiles.  It  has  been  shown  about  a  decade  ago  that 
a  computer  generated  phase  hologram,  labelled  a  kinoform,  is  entirely  equivalent 

to  a  phase  Fresnel  lens.  A  good  bibliography  on  this  subject  is  found  in 
Clair  and  Abitbol^^. 
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c)  Holographic  generation.  Instead  of  producing  a  relief  pattern  one  can 
induce  a  change  in  the  refractive  index  of  certain  materials  by  exposure  to 
light  and  realize  a  Fresnel  zone  plate.  Opterating  near  the  Bragg  angle  such 
a  Fresnel  plate  has  high  dispersion  and  focusing  power  just  like  the  BFZP. 
Recent  advances  in  the  fabrication  of  dichromated  gelatin  have  resulted  in 
very  large  index  changes  in  the  material  under  exposure  to  light.  This 
makes  dichromated  gelatin  an  attractive  material  for  infrared  applications. 
Some  problems  still  remain  to  be  solved  because  the  material  must  be  exposed 
with  visible  light  (for  maximum  sensitivity)  and  used  with  Infrared  light  for 
fiber  optic  systems.  Such  a  wavelength  shift  between  construction  and  r  con¬ 
struction  introduces  severe  aberrations. 

CONCLUSIONS 

An  off-axis  BFZP  constitutes  an  attractive  candidate  for  a  multichannel 
near  infrared  fiber  optic  demultiplexer.  More  research  must  be  devoted  to 
devise  methods  of  producing  a  BFZP. 
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SHIFTING  ZONES  (tr) 


Fig.  2  Phase  reversal  Fresnel  zone  plate  d  *  gy^y,  n  is  the 
refractive  index. 
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2-channel  8FZP  demultiplexer  with  detectors. 
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ABSTRACT 


We  formulate  a  theory  of  backscattering  from  static  perturbations  in 
fibers,  and  apply  the  results  to  predict  the  OTDR  signal  for  perturbations 
such  as  bends.  A  variety  of  effects  are  examined  in  detail,  such  as  transi¬ 
tion  loss,  curvature  loss,  dependence  on  mode  number  and  dependence  on 
pulsewidth. 
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I.  INTRODUCTION 


In  optical  time-domain  reflectometry  (OTDR),  optical  pulses  are  launched 
into  a  fiber  waveguide  and  the  back-scattered  signal  is  measured  as  a  function 
of  time  *  .  For  a  uniform  unperturbed  fiber  the  OTDR  signal  consists  entirely 
of  Rayleigh  back-scattering,  and  displays  an  exponential  decrease  as  a 
function  of  time  due  to  the  absorptive  and  scattering  losses  of  the  fiber. 

In  this  work  we  examine  backscattering  from  static  perturbations  (slowly 
varying  compared  to  the  pulse  width)  in  optical  waveguides.  The  exponentially 
decreasing  OTDR  intensity  will  in  this  case  be  modified  due  to  reflections  and 
radiation  losses  induced  by  the  perturbation.  Corresponding  changes  will  also 
occur  of  course  in  the  spatial  and  temporal  makeup  of  the  pulse  although  we  do 
not  consider  these  here. 

In  order  to  calculate  the  OTDR  return,  we  employ  a  time-dependent 
Green's  function  (GF)  method.  Although  the  equations  derived  in  this  way 
are  intractable  in  general,  it  is  possible  to  obtain  explicit  but 
approximate  expressions  for  the  case  of  weakly  perturbed  fibers  by  approxi¬ 
mating  the  exact  GF's  by  their  unperturbed  counterparts.  One  fortunate  conse¬ 
quence  of  the  lowest  order  theory  is  that  the  contributions  to  the  OTDR 
intensity  from  random  scattering  processes  (i.e.,  Rayleigh  scattering)  and 
from  deterministic  perturbations  do  not  interfere. 

The  detailed  properties  of  the  backscattered  signal  depend,  of  course, 
on  the  specific  form  of  the  perturbation  and  the  characteristics  of  the  wave¬ 
guide.  Nevertheless,  certain  general  features  emerge  from  the  present  formula¬ 
tion  as  illustrated  by  the  specific  example  of  bend  perturbations.  For  example, 
when  a  uniformly  perturbed  region  is  involved,  then  the  net  scattering  from 
the  interior  of  the  perturbed  region  vanishes,  leaving  contributions  from  just 
the  ends.  In  higher  order,  modes  may  mix  in  the  interior  and  tunnel  cut  of 
the  guide  as  well.  In  general  both  reflections  and  transition  losses  will  be 
manifested  in  the  OTDR  signal,  with  the  characteristics  of  this  signal  de¬ 
pending  on  the  symmetry  of  the  perturbation.  For  example,  for  a  single  mode  wave¬ 
guide  reflection  will  be  totally  absent  in  lowest  order  if  the  perturbation 
possesses  odd  symmetry,  but  is  present  for  even  symmetry.  Transition  losses, 
on  the  other  hand,, will  occur  in  either  instance.  We  perform  a  variety  of 
computations  to  illustrate  the  dependence  of  OTDR  from  bends  on  launch 
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conditions,  pulsewidth,  bend  length,  and  bend  radius. 

II.  GREETS  FUNCTION  THEORY  FOR  PROPAGATION  IN  PERTURBED  WAVEGUIDES 

In  this  section  we  utilize  time-dependent  GF  theory  to  obtain  an 
expression  for  the  change  in  the  time-dependent  field  amplitude  in  a  fiber 
induced  by  a  perturbation.  The  results  will  be  used  in  subsequent  sections 
to  evaluate  the  OTDR  return  from  fiber  bends. 

For  the  case  of  fiber  waveguides  for  which  the  difference  between 
core  and  cladding  refractive  indices  is  small,  which  includes  the  majority 
of  waveguides  of  interest,  the  vector  wave  equation  for  the  waveguide  modes 
may  be  approximated  by  equivalent  scalar  equations  for  TE  and  TM  modes.  Taking 
z  as  the  propagation  direction  and  assuming  solutions  of  the  form 

E(x,y,z,t)  =  E{x,y)  (1) 


the  wave  equation  for  TE  modes  may  be  written  as 


2  2  2  ) 

+  “V  +  (  n2(x,y,z)-32)/  E  (x,y)  = 

9x  9y  c  ;  y 


3y 


(2) 


For  simplicity  we  first  consider  the  case  of  planar  waveguides.  In  this 
instance  we  have 

n2(r)  =  n2(x,z)  =  nQ2(x)  +  6n2(x,z)  (3) 

2  2 
where  nQ  (x)  represents  the  unperturbed  guide  and  <Sn  (x,z)  the  perturbation. 

Eigensolutions  for  the  unperturbed  guide 
2  2  2  ) 


^7  +  (\"q2(x)  -  Snc)  \  un(x)  =  0 


(4) 


can  be  used  to  construct  the  GF  for  the  unperturbed  guide  as  c 
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-0/^  •**  \ 
G  (r,r;to) 


j 

2tt  n  )  ^2  -  a2 


dA 


+  2^df?kJ 


uk*(x/)uk(x)e'lA^z‘z  ^ 


*V 


T 


dA 


(5) 


which  can,  by  contour  integration,  be  reduced  to 


G°(r,r';u))  =  j  2  I  un*(x')un(x)e-lBJz'zl  (6) 

n  ^n 

where  the  summation  symbol  will  now  be  used  to  indicate  an  integral  over  all 
continuum  states  as  well  as  summation  over  bound  waveguide  modes.  By  con¬ 
struction  G°(r,  r';w)  satisfies  the  differential  equation 


32  I 


G°(r,  r';u>)  =  6(x-x')  &[z-z') 


(7) 


The  full  GF  G(rVr*;u) 


that  satisfies 


.2  .2  2  9 
a  .a  .  w  _  2 


2*  +  — 2"  +  n0  +  ^  5n  (x,z) 


dz 


G(r,  r^w)  =  <5(x-x')  5(z-z^) 


may  then  be  determined  by  iteration,  with  the  initial  approximation  to 
G(r,  r';<jd)  being  given  by  G°(r,  rx;w);  thus 


Gfr.r^.u))  =  G°(r,r ^oj)  +  /G°(r  r^  ;u>)  ^  5n2(r")  G°(r'i’’  rx)dr^ 

+  // G°(r?')  45n2(r^)G0(r"r^  )4>n2tf'>(^'r')dr'' dr'''  +... 
c^  cz 


(8) 


(9) 


The  scattered  wave  produced  by  the  perturbed  waveguide  in  response  to  an 
incident  wave  E°  can  then  be  expressed  using  the  full  GF  as 
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(10) 


6E(r,a)}  =  -  /  G(r  r'; u)  ^5-  <5n2(r')  E°(r/)  dr'  . 

y  cd  y 

In  order  to  exhibit  a  calculation  of  an  OTDR  backscattered  return  using 

the  GF  approach  in  the  simplest  manner,  we  shall  limit  the  present  calculation 

2 

to  planar  guides  with  perturbations  6n  (x,z)  that  can  be  factored  as 

<$n2(x,z)  =  <$n2(x)  *  f ( z )  _  (11) 


We  shall  first  assume  that  the  input  pulse  is  Gaussian  and  that  only 
one  incident  mode  i  is  excited 


.  +  i<*>  t  -  (t-t J^/zr1  -i 3 -z 

Ey  (r,t)  =  e  ui (x)e  * 

where  N  is  a  normalization  factor  for  the  input  pulse, 
normalized  to  have  equal  energy,  e.g. 

/Ey°(x,z,t)  Ey°(x,z,t)  dx  dt  =1 


(12) 


If  input  pulses  are 


(13) 


the  normalization  factor  will  be  given  by 
N  =  x-1/2  it*1'4. 


(14) 


To  calculate  the  scattered  pulse  we  decompose  the  initial  pulse  into 
its  Fourier  components  (taking  tQ  =  0) 


E  °(r,u))  =  ^  fe~,WL  E  v(r,t)dt  =  e 


1  f_-iwt  r  Of 


2  2, 


-(w-w  )  T  /2 


-  i  3  ,•  z 


/  27T 


fu^x) 


(15) 


and  propagate  each  frequency  component  ui  with  the  GF  G(r  r^m).  Using 
equation(10)and  approximating  G(r  r';w)  to  lowest  order  by  G°(r,  r*;u>),  the 
scattered  pulse  can  be  expressed  as 
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(16) 


.  .  2  .  -i6jz-z 

6Ey(r,t)  =  ///  e1wt  du>  r.  gL  ur*(x  )  uR(x)  e 


2ic  n  pn 


2  ,  -(w-w  )2t^2  ,  ^ 

•5n^(x)f(z  )e  xu.;(x  )  e  Nr  dx  dzf 

- m -  1  r 

Expanding  each  6n(w)  around  the  central  frequency  u>  as 
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3n(w)  “  6n(wo)  +  3^0,  *  (aj-wo)  +* 

o 


(17) 


=  8„°  +  v  ^  •  (lO-ld  )  +... 

n  n  o 


and  retaining  only  8n(<D0)  and  the  first  derivative  term 
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n 


3u) 


(18) 


where  vR  is  the  modal  group  velocity,  the  expression  for  the  scattered  wave 
6Fy  can  be  written  as 


Z 

n 


iw0t  i  iS^  z  .  i(u>-w0)t 
e  je 


-(w-w  )  t2/2  T  u 

a  _ — 5~ 

72tt  2ic 


+  ivn_1(o)-u)0)z 
e 


en° 


'in 


1<Bi  iBn>  Nt 


(19) 


where  e-n  and  n  are  the  transverse  and  longitudinal  matrix  elements  of  the 
perturbation  given  by 


e 


in 


=  / 


u^(x)  dx 


(20) 
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and 


<*,  +  3  )  z 

n(6ii3n)=/e  1  n  f(z)dz  (21) 

—  oo 

and  where  refer  to  backward  and  forward  scattered  waves  respectively. 

The  perturbations  considered  here  lead  in  general  to  reflection, 
radiation  loss,  and  mode-mixing.  To  obtain  predictions  of  the  backscattered 
OTDR  signal,  we  must  also  account  for  several  other  sources  of  attenuation 
in  the  fiber  such  as  Rayleigh  scattering,  intrinsic  and  extrinsic  absorption, 
and  curvature  radiation  losses  (to  be  discussed  in  section  Illb).  All  of  the 
latter  may  be  incorporated  into  the  GF  approach  by  allowing  the  propagation 
constant  for  each  mode  to  have  an  imaginary  part 


proportional  to  the  total  attenuation 


an(z)  "  an(z) RAYLEIGH  +  ^i^ABSORB  *  “n^BEND'  ^ 

The  Rayleigh  and  absorption  losses  are  nearly  uniform  throughout  a  given 
fiber  while  the  curvature  radiation  loss,  when  applicable,  is  generally 
confined  to  a  limited  length  of  fiber. 

lU.  MODEL  CALCULATION  OF  OTDR  FROM  BENDS 

In  this  section  we  apply  r.he  formalism  of  section  II  to  predict  the 
OTDR  return  from  bends  in  fibers.  For  purposes  of  illustration  we  will  simplify 
the  calculation  by  utilizing  as  our  model  a  step  index  planar  guide  with  a 
perturbation  consisting  of  a  bent  region  (possibly  containing  multiple  loops) 
extending  from  z  =  a  to  z  =  a+L  with  a  bend  radius  of  curvature  RQ.  Use  of 
the  planar  geometry  allows  us  to  demonstrate  the  principal  physical  effects 
associated  with  OTDR  from  bent  fibers  with  a  minimum  of  mathematical 
complexity.  The  methods  employed  are  easily  generalized  to  circular  fibers 
at  the  expense  of  more  tedious  calculation. 
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For  definiteness  we  choose  the  following  values  of  parameters  for 
the  unperturbed  guide:  refractive  index  of  the  core  =  1.51;  refractive 
index  of  the  clad  (assumed  to  be  of  infinite  extent)  =  1.50;  core  width 
d  =  10  pm,  center  wavelength  for  the  input  pulse  =  l^n.  With  this  choice 
of  parameters  there  will  be  a  total  of  seven  bound  modes. 

For  gentle  bends  (R  »  d)  a  bent  section  of  guide  may  be  considered 

^  O 

as  an  equivalent  straight  guide4  with  a  perturbation  given  by 


5n^(x)  =  2  x  nQ2(x) 

f0(z)  =  yiT  =  i^e(2-a>  e(«-L-z). 

For  this  case  eqp.  (19)  yields 
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(where  we  have  omitted  writing  the  superscript"o"on  3® * s ) .  From  this 
expression  it  can  be  seen  that  backward  or  forward  scattered  waves  will  appear 
to  come  from  the  front  and  rear  of  the  perturbed  region  with  appropriate  time 
delays  and  phases. 

The  reason  no  scattered  waves  come  from  the  interior  of  the  perturbed 
region  is  that  contributions  to  the  scattered  wave  from  the  interior  cancel 
out  exactly  leaving  only  the  uncancelled  front  and  rear  contributions.  In 
general  scattered  waves  will  be  generated  from  within  a  perturbed  region  if 
the  perturbation  is  not  entirely  uniform.  On  a  microscopic  scale  such  non¬ 
uniformity  is  essentially  the  mechanism  which  produce  Rayleigh  scattering. 
Inhomogeneities  in  this  case  are  due  to  frozen  in  thermal  fluctuations  in 


density  and  composition  which  result  in  fluctuations  of  the  refractive  index. 
A  calculation  of  Rayleigh  scattering  can  be  carried  out  using  the  GF  methods 
discussed  above  but  will  not  be  undertaken  here. 

When  all  processes  have  been  included,  the  OTDR  return  from  the  bend 
will  consist  of  backscattered  Rayleigh  radiation  that  has  been  diminished  by 
curvature  radiation  and  by  localized  radiation  losses  as  the  pulse  transits 
the  front  and  rear  of  the  bend  plus  a  "reflected  pulse"  that  consists  of 
backscattering  into  bound  modes  that  will  also  appear  to  come  from  the  front 
and  rear  of  the  bend.  Due  to  the  random  nature  of  the  frozen  in  fluctuations 
giving  rise  to  Rayl eight  scattering,  the  reflected  pulse  and  the  Rayleigh 
backscattered  pulse  will  not  interfere.  This  may  be  seen  by  expressing  the 
total  backscattered  power  as 
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Since  the  thermal  average  of  <5E  is  zero  .  contributions  to  the  OTDR 

n  FLT 

backscattered  pulse  from  Rayleigh  scattering  and  from  the  reflected  pulse 
will  be  simply  additive 

PT"0T(t)  =P‘R>»  +PUt>  {m 

In  the  following  sections  we  will  discuss  each  of  these  contributions 
separately.  In  section  IV  numerical  results  using  the  above  model  for  various 
curvature  radii  and  pulse  lengths  will  be  discussed. 


54 


Ilia.  TRANSITION  LOSS 

The  relative  power  that  is  radiated  away  as  the  input  pulse  passes 
the  perturbed  region  can  be  calculated  as 


where 
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(29) 


4 

are  the  forward  and  backward  radiation  losses  and  6Ey  (k)  is  the  integrand 

of  eqn.  £5)  for  the  continuum  contributions.  Due  to  the  denominator 

(0.  ±0,)~^  in  eqn.  £5)  backward  radiation  will  be  much  smaller  than  forward 

•  K  .j. 

radiation  and  can  be  neglected.  Using  the  expression  for  (k)  for  forward 
scattering  from  eqn.(25)and  performing  the  x  integration  then  yields 
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and  g. 

represent  the  radiation  pulses  coming  from  the  front  and  rear  of  the  bend. 


(31) 
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If  we  consider  the  L  dependence  of  the  forward  radiation  loss,  there 
are  two  limiting  regimes.  In  the  first  L  »  (B..-^)"1  for  all  6k  in  the 

range  n.,ko  >  0^  >o.  For  this  case  the  rapidly  oscillating  cross  term  pro¬ 
proportional  to  cos ( ( 0-j-3k)l)  can  be  set  equal  to  zero  and  the  forward 

radiation  can  be  considered  as  being  lost  in  two  steps,  the  first  as  the 
pulse  transits  the  beginning  of  the  bend  at  z=a  and  the  second  as  it  passes 
the  end  of  the  bend  at  z=a+L.  For  this  case  (after  factoring  out  the  bracketed 
term  representing  the  pulse  shape)  the  transition  loss  T.  is  given  by 
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w  1  p  3 

In  the  second  regime  L(B^-3k)~  £10  -10  the  cross  term  may  not  be 
neglected  and  there  is  the  possibility  of  radiation  fields  from  the  front 
and  back  of  the  bend  interfering.  For  this  case  f-k<u  9ik  and  the  integration 
over  all  radiation  modes  is  given  by 
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An  approximate  value  for  this  integral  may  be  obtained  by  assuming  that  the 
slowly  varying  second  factor  (hF)  can  be  taken  to  be  roughly  constant  over 
the  region  where  the  first  factor  is  a  maximum  (3-  3*»>n2kg).  With  this 
assumption  we  obtain  an  upper  bound  for  the  transition  loss 

k  4  r 

2Ti  TTY  l^6!-972  -A(1“c2(6iL)-S2(6iL)) 


-B(C2(6.L)-S2(6.L))  -  Ce  (sin^.L)*  cos^L))  + 


For  the  waveguide  we  have  chosen  as  an  example  the  6.'s  are  of  the  order 

4  -1  1 

<5.j  'v  10  m  .  The  dependence  of  on  mode  number  will  thus  be  given 

roughly  by  6^  /2. 

The  calculation  of  transition  loss  given  above  is  idealized  to  the 
extent  that  we  have  assumed  an  infinite  cladding.  For  a  finite  clad  any 
radiation  that  leaves  the  core  region  is  effectively  absorbed  locally.  The 
transition  loss  formula  without  the  interferences  term  given  by  eqn.  (25) 
would  thus  be  more  appropriate.  Likewise  the  assumption  of  the  bent  guide 
as  an  equivalent  straight  guide  obscures  the  fact  that  a  guide  may  satisfy 
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the  condition  for  interference  to  be  important  while  the  front  of  the  bend 
may  be  topologically  "far"  from  the  rear  of  the  bend  due  to  looping.  An 
additional  constraint  that  L  «  R  should  thus  be  imposed  to  further  limit 

0  A 

the  domain  where  the  interference  term  is  important. 

Illb  OTDR  FROM  BENDS  INCLUDING  CURVATURE  LOSS 

Before  the  Rayleigh  scattering  contribution  to  the  OTDR  return  can  be  obtained 
for  bend  perturbations,  we  must  first  incorporate  an  additional  loss  mechanism  due 
to  curvature  radiation  (tunneling  loss)  that  is  not  accounted  for  within 
lowest  order  GF  theory  and  that  must  be  specifically  included  in  the 
imaginary  part  of  each  propagation  constant  8n.  The  origin  of  this  curva¬ 
ture  loss  is  discussed  below. 

If  we  view  the  refractive  index  profile  as  a  quantum  mechanical  poten¬ 
tial  =  -ko2n2(x)  then  the  perturbation  <5V2  *  -kQ2  nQ^(x)  will  transform 

the  step  index  potential  as  indicated  in  fig.  1. 

It  can  be  seen  that  all  of  the  former  bound  states  now  have  the  possibility 
of  tunneling  through  the  lowered  potential  carrier  to  continuum  states  and 
calculation  of  this  tunneling  rate  will  yield  a  curvature  loss  a1RAD  propor¬ 
tional  to  ?  ? 

-  y  Ry.-  /  B  p  ?  ?  ? 

e  where  y.c  =  ^  -n2\c  and  8  =  n^.  (36) 
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Including  both  the  curvature  loss  and  the  transition  loss,  the  total 
backscattered  signal  due  to  Rayleigh  scattering  can  be  calculated  for  one 
incident  mode  i  as 


p(-)  00  -(t-(z'-z)/v.-z7v  r/  2 
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where  T.  and  R.  are  the  transition  loss  and  reflected 
power  for  the  ith  mode,  where 


(37) 


“TOT  =  “ABSORB  +  “RAYLEIGH 

are  taken  to  be  independent  of  mode  number, 

n  4 

and  a^D  is  given  for  the  step  index  guide  as 
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(39) 


(40) 


a  geometrical  factor  to  approximate  the  fraction  of  Rayleigh  scattered  power 

5 

that  is  backscattered  into  the  guide  .  The  integrals  appearing  in  this 
expression  can  all  be  expressed  analytically  in  terms  of  error  functions. 
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From  the  expression  (39)  for  aj^  it  can  be  seen  that  the  curvature 

loss  depends  exponentially  upon  both  RQ  and  the  mode  number  n.  The  step 

index  guide  we  have  chosen  as  an  example  yields  (y^)  /\j(yQ)  /12;  tunneling 

losses  for  the  weakly  bound  modes  can  thus  be  many  orders  of  magnitude 
greater  than  those  of  tightly  bound  modes. 

I  lie  DIRECT  REFLECTIONS  AND  EFFECTS  DUE  TO  GRADUAL  BEND  ONSET 

Due  to  mechanical  stiffness,  the  transition  from  a  straight  guide  to 
a  curved  section  of  guide  will  not,  in  general,  occur  abruptly,  but  will 
change  gradually  over  some  finite  distance  6  =  2Vq  as  indicated  in  fig.  2. 

In  this  subsection  we  consider  the  effect  of  a  continuous  transition  region 
(curvature  radius  changes  from  R"^=0  to  R~^  0  over  a  finite  distance) 
on  the  OTDR  return. 


We  first  consider  a  transition  model  with 

R_1(z)  =  j  £th(q(z-a))+th(q(a+L-z))J  (41) 

which  yields  for  the  longitudinal  matrix  element 
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where  =  3^  ±  3  .  In  the  point  limit  6+0,  n-*-n0, where 
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is  the  longitudinal  matrix  element  used  in  eqn.  (25)  For  the  case  of  backward 
scattering  where  q  «  p+,  eqn.  (42)  yields 
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which  goes  to  zero  exponentially,  suppressing  reflected  waves.  Only  if 
6  <  \Q  will  backward  scattered  waves  be  significant.  For  forward  scattering 
it  is  possible  for  q  to  be  comparable  with  p_  (for  the  guide  of  section  III, 
q  =  p  a/104  =>  <$*»2n  x  100  gm)  in  which  casen  can  be  of  the  same  order  as  n. 


In  a  second  model  for  the  transition  we  take 
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which  has  the  same  q-+»  point  limit  as  the  first  model.  For  backward  scattering 
and  q  <<  p+,  the  overall  magnitude  of  n(p+)  is  down  by  a  factor  q  (?c  over 
n0(P+)»  yielding  a  drastically  reduced  reflected  pulse.  For  the  reflected 
pulse  to  be  of  the  same  order  as  the  Rayleigh  backscattered  pulse,  <s  would 
have  to  be  of  the  same  order  as  X  ,  i.e.  smaller  than  the  width  of  the  guide. 
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If  the  second  model 
soidal  bends  ie.. 


is  changed  to  allow  a  uniform  series  of  a  sinu- 


1 _ =  1_  (1-cos  q(z-a) ) 

W)  "  R0  2 


a  <  z  <  a 


N  =  1 ,  2 ,  3 _ 


we  find 
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which  for  N=1  agrees  with  the  second  model  when  L  =  2rr/q  and  is  similar  in 
form  to  the  single  bend  of  the  second  model  for  N?1 . 


IV  COMPUTED  RESULTS 

To  illustrate  the  effects  of  bend  radius  of  curvature  and  pulse  length 
on  OTDR  returns  as  calculated  for  the  model  of  section  III,  we  have  plotted 
in  figs.  3  through  10  the  OTDR  return  pulse  vs.  time  for  bends  having  curva¬ 
ture  radii  of  10  cm  and  .56  cm, and  pulse  widths  of  10  ns  and  100  ns. 

For  these  curves  we  have  taken  the  longitudinal  matrix  element  to  be  n  =  nQ 
and  have  assumt  ^  reflections  to  be  zero.  The  bend  length  L  and  the  position 
of  the  bend  onset  a  have  been  taken  to  be  L  =  3.5m  and  a  =  10  m  respectively. 
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Tile  scattering  coef fieient  for  riayleiyn  scatteriiiy  and  tiie  absorption 


coefficient  for  impurity  absorption  have  been  taken  as  ^^YLE IGH  =  ctb/|cni 
and  ^bjorb  =  db/km  respectively.  For  figures  3  through  6  the  waveguide 
is  assumed  to  be  "uniformly"  illuminated  with  an  input  field  given  by 


in 


(uo(x)+u2(k)+u4(x)+u6(x)) 


with 


a.j  =  1/2  ,  i  even 
a.j  =  0,  i  odd  . 


(49) 


The  inclusion  of  more  than  one  input  mode  in  eqn.  (4Q)leads  to  a 
modification  of  eqn.  (28  )  which  represented  the  transition  loss  for  one  input 
mode,  that  includes  new  cross  terms 
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For  the  off-diagonal  terms  in  eq.  (50)  to  be  important  the  wave  packets  for 
each  pair  of  modes,  as  given  by  eqn.  (31 )  must  overlap.  If  the  onset  length 
a  is  large  and  the  group  velocities  for  the  various  modes  are  different 
enough,  the  overlap  products  will  be  very  small.  In  order  to  keep  the  present 
calculation  as  simple  as  possible  we  have  assumed  that  the  length  a  and  the 
dispersion  in  group  velocities  are  large  enough  to  enable  the  off-diagonal 
terms  be  neglected. 

Figures  3  and  4  indicate  the  behavior  of  the  OTDR  return  for  the  case 

Rq  =  10  cm  where  transition  loss  is  the  dominant  loss  mechanism.  In  figure 

3  the  pulse  width  =  10  ns  is  smaller  than  the  temporal  extent  of  the 
2L 

bend  tb  =  —  =  35  ns  while  in  figure  4  the  pulse  width  is  larger  (we  note 
o 

that  since  the  abscissa  represents  the  round  trip  time,  the  front  and  rear 
of  the  bent  section  appear  at  t  =  100  ns  and  t  =  135  ns). 

"aeffz 

Since  the  return  pulse  behaves  roughly  as  P  «  e  we  have  also 

plotted  the  derivative  of  this  curve  to  indicate  an  effective  loss  rate 

The  curve  for  in  figure  3  illustrates  the  transition  loss  from 

the  front  and  rear  of  the  bend  for  t  ^t^. 

Figures  5  and  6  illustrate  the  case  of  uniform  illumination  for 

R  =  .56  cm  where  curvature  radiation  is  the  dominant  loss  mechanism.  Since 
o 

the  input  power  is  equally  distributed  over  the  four  symmetric  modes  while 
the  outermost  modes  4  and  6  can  radiate  their  power  away  over  a  very  short 
distance,  approximately  half  of  the  input  power  (<*>3db)  is  lost  in  a  very 
short  time  after  passing  the  onset  of  the  bend  if  the  pulse  width  is  small 
(fig.  5,  x=  10  ns)  or  over  a  proportionately  longer  time  for  longer 
pulses  (fig.  6,  t=  100  ns).  (We  note  that  the  extra  wiggles  in  the  curves 
for  u  ^  in  figures  5  and  6  are  an  artifact  due  to  numerical  differentiation 
and  would  be  absent  in  a  more  exact  evaluation). 


To  illustrate  the  dependence  of  transition  loss  and  curvature  radiation 
loss  on  mode  number  we  have  plotted  in  figures  7  through  10  the  OTDR  return 
from  the  model  guide  of  section  III  excited  by  a  pulse  consisting  of  equal 
weights  of 


uQ  and  U£  (curve  A) 

u2  and  u^  (curve  B) 

u4  and  Ug  (curve  C). 

For  the  case  of  Rq  =  10  cm,  shown  in  figures  7  and  8  where  transition 
loss  dominates,  it  can  be  seen  that  only  the  least  tightly  bound  mode  6 
encounters  an  appreciable  transition  loss  at  z=a  and  z=a+L.  The  case  RQ  =  .56  cm, 
where  curvature  radiation  loss  dominates,  is  shown  in  figures  9  ana  10.  Since 

A 

ag  ~  2.3  x  103  db/km,  RA0  ~  1.4  x  10  db/km,  and  ~  .2  db/km,onl.y 

the  two  outermost  modes  are  affected  to  an  appreciable  extent.  Mode  6  is  lost 
entirely  within  a  very  short  distance  (A>l0"^m)  while  mode  4  is  only  fractionally 
depleted  in  traversing  the  bend.  The  sharp  dropoff  in  curve  C  of  figs.  9  and  10  at 

t  =  100  ns  is  due  to  the  rapid  loss  of  mode  6  input  power. 

In  figures  11  and  12  we  have  illustrated  the  OTDR  return  for  the  case 
where  n  =  nQ  and  reflected  pulses  are  not  suppressed. 

V  CONCLUDING  REMARKS 

We  have  presented  a  theory  of  OTDR  from  weak  static  perturbations  in 
optical  fibers,  utilizing  a  time-dependent  Green's  function  rormulation.  The 
formulae  we  have  obtained  incorporate  the  effects  of  direct  reflection  and 
transition  loss.  The  theory  also  incorporates  uniform  attenuation  from  absorption 
and  scattering,  as  well  as  curvature  radiation  losses  (when  applicable).  We  have 
applied  the  theory  to  examine  the  OTDR  signal  (and  its  time  derivative)  from 
fiber  bends  as  a  function  of  bend  length  and  radius,  and  the  dependence  on  mode 
number  and  pulsewidth. 
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Our  calculations  for  a  simplified  slab  guide  model  demonstrates  that 
transition  losses  dominate  for  large  radius  bends,  while  curvature  losses 
dominate  for  small  radius  bends.  In  either  case,  if  the  spatial  pulse  width  is  much 
larger  than  the  length  of  the  bent  region,  the  loss  shows  up  as  a  single  drop 
in  the  OTDR  signal.  As  expected,  the  size  of  the  loss  is  very  sensitive  to 
mode  number.  When  the  pulse  width  is  short  compared  to  the  bent  region,  separate 
features  due  to  the  front  and  end  of  the  bend  may  be  observed  in  certain  instances. 
Moreover,  if  the  bend  region  occurs  abruptly,  then  peaks  due  to  reflections  can 
be  observed  in  OTDR  with  short  pulses.  If  the  bend  onsets  gradually,  however, 
then  reflection  effects  will  be  suppressed. 

Experiments  to  date  have  indicated  agreement  v/ith  several  general  features 
mentioned  above,  such  as  the  overall  loss  at  the  bend  and  the  sensitivity  on 
the  spatial  mode  patterned  launched  in  the  fiber.  However,  further  refinements 
in  both  the  theoretical  model  (extension  to  circular  fibers)  as  well  as  the 
experimental  techniques  (better  signal-to-noise  ratio,  availability  of  shorter 
pulses)  are  required  before  a  quantitative  comparison  between  theory  and 
experiment  can  be  carried  out. 
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RELATIVE  BACKSCATTERED  POWER  (db) 


3.  Upper  curve  is  the  Log^  of  the  intensity  of  the  OTDR  return 

pulse  versus  time  for  "uniform  excitation";  lower  curve  is  the  der- 

•1  -1  3 

ivative  of  the  upper  curve  in  units  of  m  (m  *  4.343  x  10  db/km). 
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Fig.  5. 
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Same  as  fig.  3,  but  for  other  choices  of  parameters. 
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RELATIVE  BACKSCATTERED  POWER  (db) 
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Log^g  of  the  intensity  of  the  OTDR  return  pulse  vs  time:  Curve  A 
results  from  initial  excitation  of  modes  0  and  2  (in  units  of 
.016  db);  Curve  B,  modes  2  and  4  (in  units  of  .016  db);  curve  C, 
modes  4  and  6  (in  units  of  .016  db). 
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Fig.  8. 


Same  as  fig.  7,  but  for  other  choices  of  parameters. 
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Fig.  9.  Same  as  fig.  7,  but  for  other  choices  of  parameters.  Also  units 
are  changed  as  follows: 

curve  A  in  units  of  .016  db,  curve  B  in  units  of  1.6  db, 
curve  C  in  units  of  20.  db. 


RELATIVE  BACKSCATTERED  PWR  (db) 


Fig.  10.  Same  as  fig.  7,  but  for  other  choices  of  parameters.  Also  units  are 
changed  as  follows: 

curve  A  in  units  of  .016  db,  curve  B  in  units  of  1.6  db, 
curve  C  in  units  of  3.2  db. 


!  I 


Fig.  ll.  Same  as  fig.  3,  but  with  direct  reflection  included  (as  well 
as  other  choices  of  parameters). 


8 


RELATIVE  BACKSCATTERED  POWER  (db) 


R=  1.7  cm 
r=  10ns 


Fig.  12.  Same  as  fig.  7  but  with  direct  reflection  included  (as  well  as 
other  choices  of  parameters).  Also  units  are  changed  as  follows 
curves  A  and  B  in  units  of  .016  db,  curve  C  in  units  of  1.6  db. 


9 


\  I 

5  MISSION  5 

§  Rome  Air  Development  Center  ^ 

g  RA VC  plans  and  executes  research,  development,  tut  and  K 

*»  selected  acquisition  programs  in  support  oi  Command,  Control  Q 

6  Cotmunications  and  Intelligence  (C3 1)  activities.  Technical  7 

)g  and  engineering  support  within  areas  o£  technical  competence  C 
^  is  provided  to  ESV  Program  0  Alices  ( POs )  and  other  BSD  > 

jj  elements.  The  principal  technical  mission  areas  axe  \ 

?  communications,  electromagnetic  guidance  and  control,  sua-  <$ 
5  veillance  o&  ground  and  aenospa ce  objects,  intelligence  data  \ 
\  collection  and  handling,  information  system  technology,  v 

ft  ionospheAic  propagation,  solid  state  sciences,  microwave 

Sg  physics  and  electronic  reliability,  maintainability  and  rj 

compatibility.  <8 


1  b*^i*'**^*^*^>^^ 


Printed  by 

United  State*  Air  F«m 
Hantcwn  APB,  Mat*.  01731 


